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Abstract 

The purpose of relative radiometric calibration is to eliminate the strip phenomenon on the image. A new method with all 
detectors viewing the same scene is carried out. This paper provides the modeling and computation of satellite Yaw Angular 
Velocity(YAV) based on all the CCD detectors imaging the same region. The space coordinate transform process from object 
point on the earth to image point in the focal plane is found. Furthermore, the calculation method and process of YAV for 
satellite CCD camera imaging the same scene is given. Computation results indicate that the ground trajectory interval width is 
less than one pixel near the equator, and about four pixels below 35 degrees during calibration process more than 30 seconds 
imaging for 30000 detectors imaging the same scene. 
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Introduction 

At present, the optical remote sensing satellite mainly uses the linear CCD camera to obtain the remote sensing 
images. Due to the non-uniformity noise and dark current noise presented in the CCD detector, the obvious strip 
exists in the raw images obtained from the linear CCD (Corsini,2000). Relative radiometric calibration is used to 
determine the normalization correction coefficient in order to remove strip noise. 

The arrays are slightly tilted so as to reduce the effect of distortion during operation in TDI mode. For example, the 
focal plane of Pleiades-LIR is composed of five slightly tilted TDI arrays (Greslou,2012). At present, laboratory 
integrating sphere calibration method can be satisfied with the image correction before satellite launching. 
However, the calibration coefficients of the camera will change with satellite launch, aging (Markham, 1998). The 
methods of histogram matching (Horn, 1979), moment matching (Gadallah,2000) are mainly based on the 
assumption that each sensor is exposed to scene radiance with approximately the same probability distribution. 
The side-slither maneuver method is a new method, which involves yawing the satellite 90°(Cody,2011). This 
paper provides the computation method of satellite YAV for slightly tilted CCD viewing the same scene, based on 
the establishment of strict conversion from geographic coordinate to focal plane coordinate. 

Coordinates Transformation for Satellite Squint Imaging 

There are seven coordinates involved in the coordinate transformation from object point on the earth to image 
point in the focal plane. They are defined as follows: 



Focal Plane Coordinate System (FPCS): P0-P1P2P3 

Origin Po is center of CCD arrays, and Pi, P2, P3 are parallel to roll, pitch and yaw axes, respectively. As shown in 
FIG.l, CCD camera is in conventional imaging mode, the tilted five CCD arrays are closer to being perpendicular 
to the satellite flight direction, and the incline angle with each other is 0.3°. When on-orbit running, considered 
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influence of lens distortion and earth curvature, CCD arrays should be carried with geometry calibration for 
obtaining the actual model. Here the five CCD arrays model is fitted with circle. Poins A,B,C,...,J are the end and 
mid points of five arrays analyzed in after experiments. 

Geographical Coordinate System (GCS): G0-G1G2G3 

Go is the origin of GCS, and also the corresponding ground point of Po. Gi is parallel to flight direction. G 3 is 
pointing to zenith. G 2 is perpendicular to G 1 G 3 plane. 

Inertial Coordinate System (ICS):0-I1I2I3 

Origin O is at geocenter, I2 points to north pole, I3 pointing satellite orbit plane and the equatorial plane 
intersection, Ii is perpendicular to I2 and I3. 

Earth Fixed Coordinate System (EFCS):0-E1E2E3 

E 2 points to north pole. It is in inertial space, around the counter clockwise direction with velocity 00 . 

Orbit Coordinate System (OCS):BO-B1B2B3 

Origin Bo is located on the orbit, and Bi points to the track direction, B 2 is perpendicular to the surface of the 
satellite orbit, B 3 towards the zenith. 

Body Coordinate System (BCS):S0-S1S2S3 

Origin is the same as that of OCS. The three axial positions of 4> ,0 ,i[j mean that of BCS in OCS. 

Camera Coordinate System (CCS):C0-C1C2C3 
Without installation error, CCS is the same as BCS. 

This proposed transformation model takes into accout not only 4>,0 ,i|j equal to zero, but also unequal to zero, 
which is different with Ref.(Wang,2005) only considering nadir point. Transformation from object point in GCS to 
image point in FPCS is established: 

Following the translation procedures and relationship shown in FIG.2 and FIG. 3, the transformation from GCS to 
FPCS is expressed as the following formula (1): 
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a) 


Where, R is earth radius, H is satellite height, h is object height, io is orbit incline angle, y is geocentric angle 
between satellite and ascending node, y=yo+Qt, Q is satellite angle velocity, 00 is earth rotation velocity, i[;=i[Jo+di[T, 
0=0o+d0t, cj)=cj)o+dcj)t, i[jo, 0o, cf>o are initial angles, and di[),d0,dc|) are velocity. ZK 0 OK 1 and ZK 1 OG 0 are geocentric 
angles with squint imaging mode, shown in FIG. 3 (right), which are computed: 

In ABqOKj, ZK 0 OK 1 =arcsin[(R+H)sin(p/R]-(p, B 0 K, = /(R + H) 2 +R 2 -2-(R + H)-R cosZK 0 OK, 
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( R+H I 

InAB 0 OGo, ZOB o G o =arccos(cos0coscp) , ZB 0 OG 0 = arcsin 1 sm ZOB 0 G 0 J-ZOB 0 G 0 , 

B 0 G 0 = ^/(R+H) 2 + R 2 - 2 ■ (R+H) ■ R • cos ZB 0 OG 0 

In ABoKiGo, KjG 0 = /b o Kj 2 + B 0 G 0 2 2- B 0 Kj ■ B 0 G 0 -cosO , Therefore, in AK 1 OG 0 , ZKjOG 0 = 2 -arcsin ^ 0 
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FIG. 2PROCEDURES OF TRANSFORMATION 




FIG. 3 CONVERSION FROM GCS TO FCS (LEFTiREL ATION, RIGHT: GEOCENTRIC ANGLE CAUSED BY SQUINT IMAGING) 


Computation of Satellite YAV 

As shown in FIG.4, in orbit position A, satellite is with conventional imaging mode, according to the formula (1) 
the drift angle [3o can be calculated. And then the satellite yaw angle is rotated 90- |3o degrees clockwise direction, 
and in orbit position B, satellite is with the relative radiometric calibration mode. Where the |3o is to remove 
influence of earth's rotation. And focal plane for calibration mode is shown in FIG.4(RIGHT). 

However, when on-orbit calibration imaging, due to the slightly tilted of CCD arrays, it is needed to constantly 
adjust yaw angle. To make all the detectors imaging the same region, they must be strictly through the same point 
in the natural scene, and all detectors along the same trajectory ground, so a large amount of data all the detectors 
imaging the same region can be obtained. YAV can be computed as following steps: 

1) As shown in FIG. 5, suppose the start time to of performing relative calibration is imaging zero. When the 
corresponding FPCS with origin Po is known as the original FPCS (OFPCS), and the corresponding GCS with origin 
Go is known as the original GCS (OGCS). 

2) OFPCS shifts along satellite flight direction with CCD size de, simultaneously OFPCS rotates aroung origin Po 
with certain YAV. After time ti, ti=ni*te, OFPCS will arrive at new point Po', which is the origin of new FPCS 
(NFPCS) Po'-PTPz'Ps'. 

Where te is the integral time of CCD, m is the moving pixel number of FPCS, therefore coordinate value of Po' in 
NFPCS can be computed with p 10 ^=n 1 d c sinp 0 , p 20 -n 1 d c cosp 0 . According to formula (1), the corresponding 
coordinate value of Po' in OGCS can be calculated, defined as Go'(gio',g2o',0). 

3) According to the condition of all detectors viewing the same ground scene, suppose point N in OFPCS and point 
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M' in NFPCS pointing to the same ground point, as shown in FIG.5. Set the coordinates of point N in OGCS is 
expressed as PN(piN,p2N,0), and compute corresponding geographical coordinates G N (gi N ,g 2 N,0) in OGCS according 
to formula (1). 


Suppose coordinates of point M' in NFPCS is PMinew(piM'_new,p2M'_new,0). The relationship of M' and N in FPCS is 
p 2 M'_new = (p 2 N-p 2 o V cos|3o), likewise, the corresponding coordinates of M' in NGCS is computed, expressed as 

Gm'_ _new (glM'_ .new, g2M'_new, 0) . 


4) Then , with the condition of viewing the same scene, there is g2M'_ new +g2o' equals to g2N. Get the complex 
trigonometric function equation, and the satellite YAV can be calculated. 



FIG. 4 ON-ORBIT CALIBRATION (LEFT:YAW ROTATION FOR CALIBRAITON, RIGHT: FOCAL PLANE WHEN CALIBRATION) 



FIG. 5 DIAGRAM OF SATELLITE YAV COMPUTATION FOR ALL CCD DETECTORS VIEWING THE SAME SCENE 


Experiments and Analysis 

In this paper, experimental simulation parameters are: per linear CCD length 6000 pixels with 5 pieces, focal length 
2600mm, integration time 362.59ps. Satellite orbital height H 645.338km, earth radius R 6378.137km, satellite orbit 
period 98 min, orbit incline angle 98°. 

FIG. 6 shows trajectories of 11 detectors (Po, A, B, C, ..., I, and J) in 5 CCD arrays, as shown in FIG.4 (RIGHT), 
within about 32 seconds, where, at 0° latitude, YAV adjusting period equals to 0.36s. The different color curves are 
the trajectories. As distortion makes the projection of the detector unstraight, the trajectories are also curved. The 
rectangular box area with red dashed line is common area with all 11 detectors passed through, where the 
maximum mean width within the region is 2.157m. In order to see clearly this thickness, the trajectories on the 
ground are enlarged, shown in FIG. 6(b). FIG. 6 (c) is YAV curve, which is first increased from -0.06761% to - 
0.06617% and then decreased to -0.06747%. 

Table 1 gives the mean width of trajectories with YAV adjusting period all equals to 0.36s at different latitudes and 
roll angles. It can be seen that near the equator, the width is very narrow, such as width is 2.157m at 0 Q , and 1.378m 
at 5 Q . The ground track accuracy is better than 1 pixel for 2.5m resolution satellite. With the increasing of latitude, 
the ground track width is wider and wider. And until 50 degrees, the ground width has reached 15m (6 pixels). 
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With roll angle increasing, the width near equator is the most obvious changing, such as from roll angle from 0 Q to 
40 Q , width changing from 2.157m to 22.287m. 
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(b) ZOOM ON REGION 1 



FIG.6 TRAJECTORIES OF 11 DETECTORS IN 5 ARRAYS (LATITUDE=ZERO DEGREE, YAV ADJUSTING PERIOD=0.36S) 


TABLE 1 TRAJECTORY WIDTH AT DIFFERENT LATITUDES AND ROLL ANGLES (METER) 
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Conclusion 

To view the same scene for the relative radiometric calibration of multi tilted CCD arrays, this paper has 
established the rigorous mathematics model between FPCS and GCS, where the geocentric angle is computed. 
Furthermore, it has put forward the YAV calculation method. Simulation experiments give the trajectories width, 
and also analized the width changing at different latitudes and roll angles. 
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